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ABSTRACT 

A self-consistent, thermodynamic approach is employed to derive the wave energy of a mag- 
netohydrodynamic system within the harmonic approximation and to obtain the familiar 
dispersion relation from the resulting equation of motion. The evolution of the system due 
to an external perturbation is studied by a linear response formalism, that also gives the 
energy absorbed by the magnetohydrodynamic system from the external field. The calcu- 
lated wave energy reveals the presence of positive and negative energy waves, that coalesce 
together to give rise to Kelvin - Helmholtz instability of the system. The threshold value 
of this instability changes only slightly in the presence of a small amount of viscosity, thus 
precluding the dissipative instability of the negative energy waves. The prediction of such a 
dissipative instability by earlier authors turns out to be the result of an erroneous choice of 
the viscous drag force, that violates the fundamental law of Galilean invariance. 

1. Introduction 

Propagation of waves in a medium is the manifestation of the system's response to a small 
deviation from its local thermodynamic equilibrium. The dispersion relation D{k,n) = 
describes the frequency (n) vs. wavenumber (k) relationship of such waves, and the form of 
the dispersion function D{k, n) is decided by the restoring forces and the degrees of freedom 
that the system possesses. 

The degrees of freedom of a magnetohydrodynamic system consist of its density, pressure, 
velocity and its magnetic induction. The conventional method of determining the stability 



of the system and also the nature of its wavemodes, is to work with the (MHD) equations 
connecting the above variables, thus obtaining the dispersion relation by requiring a non- 
trivial solution of the problem (eg. Alfven 1950; Cowling 1957; Chandrasekhar 1961). This 
method, though adequate to study the system's stability, does not directly allow us to 
calculate the energy of the system. 

In this paper, we study the response of a magnetohydro dynamic system to an exter- 
nal perturbation from the point of view of its energetics. The system consists of a two- 
dimensional slab-like inhomogeneity in the plasma and the flow parameters of an otherwise 
uniform magnetic medium, so that the equilibrium pressure, density, magnetic field and 
velocity are given by 



Po{z), po{z), Bo{z), uo{z) 



POo, POo) Bqo, Uqo, \z\ < d, 



POe, POe, Boe, Uoe, \z\ > d, 

with the magnetic field vectors and the steady fiows being aligned to the axis of the slab, ie., 
to the X-direction.The equilibrium condition demands that the plasma inside and outside 
the slab are in the total (gas -|- magnetic) pressure balance, namely, poo + Bq^/Stt — poe + 
Bq^/Stt. We here note, that Nakariakov and Roberts (1995; see also, Satya Narayanan 1991; 
Nakariakov, Roberts and Mann 1996) analysed the normal modes of the above equilibrium 
through the solutions of the slab's dispersion relation. These authors also included plasma 
compressibility in their analysis. Considering the algebraic complexity of the calculations of 
the slab's energetics, we here confine ourselves to an incompressible slab. 

In what follows, we begin with the basic definition of the energy density of a magneto- 
hydrodynamic system, deriving ultimately an expression for the space-averged total energy 
{Hamiltonian) of the perturbed MHD slab in terms of its generalized co-ordinates, namely, 
its interfacial displacements r]{x,z — ±d,t), and their time derivatives f]{x,z — ^d,t); see 
Section 2. In Section 3, we derive the equation of motion of the perturbed slab from the given 
Hamiltonian, which we solve directly to obtain the energy of excitation of the wavemodes. 

The above calculation of the wave energy reveals the presence of negative energy waves 
(NEW), the existence of which was predicted earlier by various authors in both the ordinary 
hydrodynamic (eg. Benjamin 1963; Cairns 1979; Craik and Adam 1979; Ezerskii, Ostrovskii 
and Stepanyants 1981; Ostrovskii and Stepanyants 1982; Craik 1985; Ostrovski, Rybak and 



Tsimring 1986) and the magnetohydrodynamic (eg. Acheson 1976; Ryutova 1988; Ruderman 
and Goossens 1995; Ruderman et al. 1996; Joarder, Nakariakov and Roberts 1997) systems. 
The present paper improves upon these earlier calculations by extending the results of Cairns 
(1979) to the much complicated magnetohydrodynamic situations, thus providing a simple 
expression for the wave energy in terms of the hnear dispersion function D{k, n) of the MHD 
slab. In his derivation of the wave energy, Cairns foUowed a procedure that is somewhat 
similar to the ones developed earlier by Stix (1962) and Witham (1974), and used the 
expansion pi_2 ~ -Di,2 {ojq + i^, /co) -^{t) exp {ik^x — out) (see. Cairns 1979) for the pressures 
across the interface of the fluid. Such an expansion, though consistent with the dispersion 
relation, appears to us to be intuitively presented by Cairns. Our method, on the other 
hand, gives rigourous calculations of all the physical parameters and, in addition, allows one 
to calculate directly the energy of the system (Section 2), from which both the dispersion 
relation and the wave energy follows, via. the equation of motion (see Section 3). 

Along with many other properties, the negative energy waves are also supposed to ex- 
ihibit dissipative instability, so that the waves become overstable in the presence of any 
arbitrarily small amount of viscous dissipation in the medium that is in the rest frame of 
these negative energy waves; cf. Kikina (1967); Weissman (1970); Cairns (1979); Ezerskii 
et al. (1981); Ostrovskii and Stepanyants (1982); Craik (1985); Ostrovskii et al. (1986); 
Ruderman and Goossens (1995); Ruderman et al. (1996). In this paper, we show that this 
viscous overstability of the NEW is simply the result of an erroneous choice of the dissipative 
damping, that violates the law of Galilean invarience. In Section 4, we show that both the 
positive and the negative energy waves, in fact, exhibit dissipative damping under the ac- 
tion of a viscous drag force, that is consistent with an appropriate Galilean transformation. 
Concluding remarks are given in Section 5, indicating the relevance of this study to several 
astrophysical MHD systems. 

2. Total energy of the perturbed slab 

2.1. The Perturbations: Consider then, that due to the action of some external stress, the 
interfaces z = ±d of the magnetic slab (see the previous section) are displaced by an amount 
r]{x,z — ±d,t) — fi{k,z — ±d,t) exp{ikx), where fj{k,z — ±d,t) are the amplitudes of 
the Fourier components of the displacements r){x, ±d, t) with respect to x, and k is the 
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wavenumber. The associated z-displacements of the media internal {\z\ < d) and extrenal 
(1^1 > d) to the slab are given by 

f}{k, z — —d, t) exp (ikx) exp {k{z + d)}, z < —d, 
r){x,z,t) = { [acosh{kz) + Psmh{kz)]cx-p{ikx), \z\ < d, (2) 

fj{k,z = d,t) exp {ikx) exp {—k{z — d)}, z > d, 
thus showing that for an incompressible plasma slab, as is the case considered here, the 
perturbation amplitudes fj{k,z,t) are evanescent in both the internal (|2;| < d, denoted by 
'o') and the external {\z\ > d, denoted by 'e') media away from the interfaces z = ±d of the 
slab. Such perturbations pertain to the surface modes of the slab with the 'cosh' solution 
in equation (2) presenting a kink (even) surface mode, and the 'sinh' solution presenting 
a sausage (odd) surface mode; cf. Roberts (1981a, b), see also Edwin and Roberts (1982, 
1983). 

The coefficients a and P in equation (2) can be evaluated from the condition of continuity 
of the z-dispacements across the interfaces z = ± d oi the slab. Retaining the above 
terminology to reperesent the perturbations that are of even or of odd symmetry with respect 
to the axis {z — 0) of the magnetic slab, we thus obtain 

a = fiy,iny,{k,t)/ cosh{kd), fii,ini,{k,t) = l{fi{k,z ^ -d,t) +fi{k,z = d,t)} , 
P = 'nscLViscLge{k,t)/ sinh {kd), r]sa.nsa.ge{k,t) = ^ {f]{k, z = -d,t) - r]{k, z = d,t)} . 
So far, we considerd only the z-displacements of the magnetic plasma. The perturbations 
in other magnetohydrodynamic variables can be obtained by applying equations (2) and (3) 
to the linearized MHD equations (eg. Chandrasekhar 1961). In medium 'o' internal to the 
slab, these fluctiations are written as 



f / 7 \ f~ /7 . sink (kz) _ cosh(kz)] .., . . 

5u^o = - (n + kuoo) |?7kink(A:, t) ^^^^ ^^^^ + ?7sausage(fe, t) ^.^^ ^^^^ | exp {tkx), (4 a) 

. fn + kuoo\ [~ ,-,cosh(kz) ,smh.(kz)] ... , 

Su^o = ^^^J|^>ci„k(^,^)^^^ + W(^'^)-^|exp(zM, (4 b) 

rr-. r i ~ /, . Smh.( kz) ^ . COSh (/cz) 1 /., x r,^ 

5B^o = -^^oo|r;.,.(fc,t)^^^ + Wge(^,^)-^|exp(zM, (4 c) 
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and 

( kBQo\ { cosh, ikz) ■ . sinh (kz)\ / , x /. ,n 

Here, ^m^: and 6uz are the perturbations in the x- and the z- components of the velocity, 
whereas, SB^ and SB^ are the perturbations in the x- and the z- components of the magnetic 
field, respectively. In writing equation (4), we have used only the n-th temporal Fourier 
mode, so that fj{k,t) = infi{k,t). For any arbitrary fj{k,t), equation (4) can be generalized 
by appropriately summing over all possible n s; such an exercise, however, does not alter the 
results that follow, as can be verified by using the condition of independence of the Fourier 
components. 

The fluctuations of the MHD variables in the external (|2;| > d, denoted by 'e' ) media 
can similarly be derived. They are 



and 



5uxe — ^{n + kuoe)f]{k, z — ^d,t) exp{ikx) exp {±k{z ± d)}, (5 a) 



Suze — ( j fj{k, z — ^d,t) exp(ikx) exp {±k{z ± d)}, (5 b) 



SB^e = T {kBoe) r]{k, z = ^d,t) exp (ikx) exp {±k{z ± d)}, (5 c) 



SBze — (- — —jr]{k,z = ^d,t)exp{ikx)exp{±k{z±d)}. (5 d) 



2.2. The Energy Densities: To calculate the total energy (or, the Hamiltonian) of the mag- 
netic slab, we first note that, the energy density of an incompressible magnetohydrodynamic 
system is given by (see the Appendix) 

where, p = po + 6p, B = |B| = |Bo + 5B|, and u = \u\ = |uo + Su\ are the total pressure, 
the total magnetic field strength and the total velocity of the medium, respectively. For the 
specific situation considered here (see Section 1), we can now use equation (6) to write the 
magnetohydrodynamic equations of motion in the following particular form: 



V£ = -ipo{n + kuo)5u + —^SB + —V\5B\'^ + PoUqV (Su^) , (7) 

47r OTT 

that must be integrated, with the help of equations (4) and (5), to obtain the expressions 
for the energy densities separately for each of the media internal ('o') and external ('e') to 
the slab. Consider first the internal ('o', \z\ < d) medium. Equation (4) then allows us to 
evaluate the various terms on the left hand side of equation (7) upto the second order of 
smallness in f}{k, t). Writing r] (x, t) — fj {k, t) exp (ikx), we thus get, 



Stt 



sinh (kz) 



x,t) 



. . cosh (kz) . , sinh {kz 

^" ^sausage 



cosh {kd) 



sinh [kd) 



cosh {kz) 
sinh {kd) 

+ 



POo" 



{n + ku. 



OoJ 



k 



An 



^sausage (-^j ^) 



cosh 



, . sinh {kz) 

POoUOo {SUxo) 



(Ml 



+ 



(8) 



sinh {kd) j 

which, after integration, gives the following expression for the energy density Eq of the 
medium internal to the slab: 



POo 



{n + kuooY 



kB, ^ 



-'Oo 



k 



An 



sinh {kz) 



cosh {kz) 1 



%ink(x,t)^^^^^^^^ +77sausage(^,t)g.^j^^^^^ | 



+ 



Stt 



sinh. {kz) . . cosh {kz) 

^kink(2^) t) 1 / , ,s ~l~ ^sausage ^)~ 



cosh {kd) 



sinh {kd) 



+ 



f}kink{x,t) 



cosh {kz) 
cosh {kd) 



'sausage 



{x^t)-: 



sinh {kz) 



+ POol^Oo {^U^o) + £0o + ^Sqo{x, Z, t) + Ko. 



(9) 



sinh {kd) 

Here, the quantities £oo and are constants of integration, that are independent of x 
and z. The function 5eoo{x, z,t) in equation (9) represents some initial fluctuations in the 
equilibrium energy density £oo — Poo + B^/Stt + (poo/2) u^^ of the medium 'o', that are the 
manifestations of some externally applied initial stress on the system. Such a stress does not 
appear explicitly in equation (7), but is incorporated phenomcnologically in equation (9) to 
facilitate the study of the system's response rj{x,z,t) to such external perturbations Seq in 
Section 3.2 of this paper. 
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In what follows, we further consider the fluctuations 6eo{x,z,t) to have only a piece- 
wise dependence on the 2;-lo cations; ie., for any particular value of x, the magnitudes of 
5eo{x,z,t) may change abruptly across the interfaces z — ±d of the slab, but are constant 
in each of the three media, namely, the medium 'o' inside {\z\ < d) the slab, and the 
medium 'e' on either {z < d and z > d) side of the slab. For the medium 'o', we can write 
Seoo{x, z, t) = Ssqo{x, z = 0,t) in such a case. 

Unlike £oo and Seoo{x,0,t), the quantity Kq in equation (9) contains terms, that are of 
second order smallness in the amplitudes of the interfacial displacements fjlx, ±d,t). Such 
second order terms do not follow from the linearized equation (8), but are ought to be 
introduced, as we have to determine all the harmonic terms in the system. A general 
expression for Kq can therefore be written as 

Ko ^ a fj^^^ik, t) + b Cusagel^^ + C fj^nkik, t)fjss.use.ge{k, t) , (10) 

where, the coefficients a, h and c must be determined separately from some physical con- 
siderations. In order to determine these coefficients, it is convenient to first drop the term 
P{k>UQo{5uxo) from equation (9). Such terms do not contribute to the total Hamiltonian 
H{fj{k,±d,t),fj{k,±d,t)) of the system, as the integrals like poUoJJ (Su^) dzdx represent 
the velocity of the centre of mass of the respective layers, and hence are identically zero 
owing to wave propagation. Unfortunately, this property of the integrals has to be imposed 
and cannot be demonstrated here, as the linearization limit that we follow, introduces some 
unphysical 77^ (/c, ±d, t) -type terms in the Hamiltonian of the system. 

In order to identify the constants a, b, c in equation (10), we use now the condition of 
equipartition. After dropping the term involving 6uxo in equation (9), we next consider an 
averaging of the energy density over a time-scale that is much longer than any periodicity 
present in the system. Such an averaging retains only the second order terms on the right 
hand side of equation (9), while all the terms linear in fj{k,±d,t) drop out. By denoting 
time averages by overbars, and by using the Parseval formula for Fourier transforms, we thus 
obtain 
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2[£o (X, Z, t) - £oo - SSoo (X, 0, t)] = 



k'B^^ sinh" (kz) 
Stt cosh^ (kd) 



^kink (^1 ^) 



+ 



Stt sinh (kd) 



+ b 



^sausage 



{k,t) + 



cosh (kz) -ji 



cosh {kdi) 



^kink 



(11) 



sinh {kz) -ji a. \ 
+ ~r~r9~rr~;7 ^sausage \ki*) 



+ 



sinh2(A;rf)''^'^"^^^^ 
A;252^sinh(2A;z)- 



+ 



sinh (2A;^) 
47r sinh {2kd) 



+ c 



sausage 



sausage 



47r sinh (2A;(i) 

On the basis of the principle of equipartition of energy in harmonic oscillators, we may now 
argue, that the contributions from the terms containing ^/^kink/sausage)(^; ^) ^lust be equal to 

. 2 

the contributions from the terms containing ?7(kink/sausage)(^' ^PPtyiiig similar arguments 
for the cross terms (containg ?7kink(^, ^)?7sausage(^) etc.), we ultimately get 



-°0o 



c = 0. 



Stt cosh2(M)' " Stt sinh^(A;rf)' " ^^^^ 

Using equations (9-12), we finally arrive at the expression for the energy density of the 
medium 'o' internal to the slab. This expression reads 



£o{x, z, t) = 6oo + Seoo{x, 0, t) + 



Stt 



sinh^ (kz) 
cosh^ (kd) 



+ 



l^?kinkOM)_ 

2 cosh^ {kd) 



+ 



.2 ^A^Osh^ (/CZ) 1 ^sausage (a^,^) 



/sausage 



[x,t) 



sinh {kd) 2 sinh {kd) 



sinh {2kz) \ k^B^ 



Oo 



+ 2?7kink(x, t)rj,,ns.se{^, t) ^.^^ ^^^^^ | + ^^J 

.2 , sinh^ (/cz) ( \- ( \ 

"I" '^sausagel-^' '^JZa^TJi^ '^^kink(^) * j ^sausage *) 



^kmk(^) ^) 



cosh^ {kz) 



+ 



POo- 



sinh {kd) 
{n + kuQof kBl^ 



cosh^ {kd) 
sinh {2kz) 
sinh {2kd)_ 



k 



in 



. sinh {kz) , , cosh (A;z) 

%ink(^) / j\ ^sausageV-^) ^j" 



(13) 



cosh {k(£) sinh (/erf) J 

The energy density of the external media can be derived by following a similar procedure. 
This energy density is given by 



Sq{x^ z, t) 



Stt 



v\x,^d, t) + —7i\x, ^d,t) 



exp{±2A; {z±d)] 



kBl 



- POe 



{n + kuQeY 



Air k 
+ eoe + Seoe{x,Td,t), 



r]{x, ^d, t){±k {z ± d)} 



(14) 
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where, the '-' and the '+' signs denote the regions z < —d and z > d, respectively. In writing 
equation (14), we have assumed that some initial fluctuations, in the form of Seoe{x, ^d,t) 
, is imposed on the equilibrium state £oe of the external plasma on both the sides of the 
magnetic slab, that are the manifestations of some external stress applied on the system; see 

the discussions following equation (9) above. 

2.3. The Total Hamiltonian: The total energy, or the Hamiltonian H[ri, fj) of the given 
MHD system is derived by integrating its energy density over the entire volume of the 
system consisting of the magnetic slab in its magnetic environment. Thus 



H (fj{k, z — ±d,t),fj{k, z — ±d,t)^ — j j e{x, z,t) dz dx 

-d+ii{k,z=—d,t) exp [ikx) pd+ii{k,z=d,t) cxp iikx) 

' ' dz 

!,t) exp (ikx) 



/r—a+ii(ii,z=—a,t)exp{iKX) na-\-ii(ii,z=a,t) cxp 

I ee{x,z,t)dz + eo(x,z,t)dz 

I ^— oo J —d+ri{k,z=—d,t) ei 

+ / eJx,z,t)dz >dx, (15) 

Jd+fj{k,z=d,t)e:x.p{ikx) J 



where, the quantities So{x,z,t) and Se{x,z,t) are as given in equations (13) and (14), re- 
spectively. 

After performing the integrations as indicated in equation (15), while retaining terms, 
that are of second order of smallness in the displacements r]{x, z, t), we finally arrive at the 
required expression for the total energy of the perturbed magnetic slab. This expression 
reads 



H (?7kink, ?7sausage, j7kink) ^sausage) — ~Crkink(— ^, ^)?7kink(^, t) - 

sausage 

{-k,t)fj 

sausage 

{k,t) 



+ %ink(^>^) 



+ 



+ 



POe 



[n + kuQ^y 



?,kBl 



^ -"Oo 

Stt 



k IGtt 
2sin\i{2kd)+A{kd) 



+ 



{n + kuooY 



Sk cosh^ (kd) 



sausage 



k 



POe 



47r 

2 



tanh {kd) 



[n + kuoe) 



3k 



POo 



[n + kuQoY 



k 



coth {kd) 4- 



Stt 



k IGtt 
2smh{2kd)-4:{kd)' 



8k sinh^ (kd) 



+ Vkink{k,t) 



167rn^ Stto^ 



+ n. 



sausage 



{k,t) 



2 sinh {2kd) + 4 {kd) 
8k cosh^ {kd) 
2 sinh {2kd) - 4 {kd) 



(16) 



S/csinh^ {kd) 
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where, we have used the definitions 



Seoe{x,-d,t) - SeQe{x,d,t) = ^^a^^i^kik' ,t) exp (ik'x), (17 a) 

k' 

and 

5eoeix, -d, t) + Seoeix, d, t) - 2feoo(2;, 0,t) = ^ 5-sausage {k', t) exp {ik'x), (17 b) 

k' 

with 5"kink {k', t) and CTgausage {k', t) being the k' -th Fourier amplitudes of the externally ap- 
plied stresses, that excite the kink (even) and the sausage (odd) perturbations of the slab, 
respectively. 

3. Excitation energy of the wavemodes. 

3.1. Equation of Motion and the Dispersion Relations: The Hamiltonian, that we derive in 
the previous section (see equation (16)), immeditaely allows us to write down the equation 
of motion for the perturbed magnetic slab. This equation reads 



2 sinh {2kd) ± 4 (kd) 



Ak cosh^ {kd) 



^(kink/ sausage 



+2 



POe 



k 



3kBl 
IGtt 



+ 



POo- 



(n + kuQoY 



k 



An 



fik) 



+ 



k'^Bl^ \2smh{2kd)±Akd 
Stt 



8A; cosh^ {kd) 



'7(kink/sausage) 



5^(kink/sausage) ( k,t) , 



(18) 



where, the '+' and the '-' sign applies to the cases of the kink (with f{k) = tanh (kd)) and 
the sausage (with f{k) = coth {kd)) perturbations, respectively. 

To check the correctness of equations (16) and (18), we consider, for the time being, 
that there is no external stress apphed on the system, ie., a{—k,t) — 0. The motion of the 
plasma slab then consists of its various eigenmodes, as revealed by its dispersion relations 
D{k,n) = 0, connecting the frequency (n) and the wavenumber (k) of any temporal Fourier 
component fj{k, n) exp (int) of the displacements fjlk, t) of the slab. For the kink and the 
sausage disturbances, these dispersion relations are obtained by substituting d/dt — in in 
equation (18). For arbitrary rj{k,t), we thus obtain 
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n + kuoo) -k'ci, f{k) = 0, (19) 



-D(kink/sausage)(/i'", n) = poe + kuoef - k^c\^ + POo 

where, we have introduced the Alfven speeds cao = -Boo/(47rpoo)^''^ and CAe = -Boe/(47rpoe)^''^ 
in the two media 'o' and 'e', respectively. 

We here note, that equation (19) is essentially the same as the dispersion relations ob- 
tained by Nakariakov and Roberts (1995), for the normal modes of an incompressible mag- 
netic slab embedded in an incompresible MHD medium, where there is a relative tangential 
velocity between the slab and its environment. This equation admits solutions 



-k 



1 /2' 



with 



a, 



e/ 



POe + POof{k) ' Poe + POof{k) ' ""'"^ "^^^ Poe + POof{k) 



Poofik) 



-,and clf = 



POcCac + Poof{k)c 



Ao 



(20 a) 



(20 b) 



Each of the kink {f{k) — tanh (/co?)) and the sausage {f{k) — coth{kd)) solutions then 
allows two distinct eigenmodes of the slab, that are represented by the and the '-' 
signs in equation (20). These modes are purely oscillatory (giving surface waves) when the 
discriminant of equation (20) is real, but one of them (the '+' mode) becomes a growing 
mode (giving Kelvin-Helmholtz (K-H) instability) when this discriminant is imaginary. The 
instability thus sets in for relative velocities 



koe - > C),f/ (ao/CVe/) ' , (21) 

ie., for values at which the '+' and the '-' surface modes of either the kink or the sausage- 
type coalesce together in the real n vs. k plane to produce an unsatble region of complex 
frequency; see Cairns (1979) for some examples of such coalescence instabilities drawn from 
hydrodynamics. Such a coalescence instability occurs only when the modes involved have 
energies of opposite sign; see the discussions in the next section. We also note, that equation 
(21) reduces to the instability criterion given by Singh and Talwar (1994) and Nakariakov 
and Roberts (1995) in a situation, where the plasma slab moves in a static environment, ie., 
when uoe — 0. 
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Equation (21) further shows, that the kink and the sausage modes merge together to give 
only two ('+' and '-') surface modes in the case of an infinitely wide {kd oo, f{k) 1) 
slab, so that the quantity Ck/ coincides with the phase speed — (poocio + Poecie)^'^^/ (Poo + Poe) 
of the hydromagnetic surface waves (Roberts 1981a,b) in a single surface of discontinuity sep- 
arating two uniform magnetic plasma media. The instability criterion reduces to the classi- 
cal threshold for the K-H instability of a magnetic tangential discontinuity (Chandrasekhar 
1961) in this limit of an infinitely thick plasma slab. 

3.2. Evolution of the System under External Stress: With this brief discussion on the nature 
of the normal modes of the slab, we now turn to the calculations for the work done by an 
external stress a{—k,t) in exciting each of these modes, which, in turn, is stored as the 
energy of that particular wavemode of the system. To find this wave energy, we consider 
that the external stresses begin to act on the system at a time t — 0, so that 



^■(kink/sausage) (-^,^) = 5-(kink/sausage) (-^)0(^)-P'(^) , (22) 

with F{t) being an arbitrary function of time t, and Q{t) being a Heaviside unit step function. 
The causal response of the system to this external stress is given in terms of a response 
function G{k,t). Thus 



^(kink/sausage) 

{k,t) 



if i < 0, 

J (kink/sausage) {k, t — t') F{t') dt' if t > 0, 

the Laplace transform of which is given by 



(23) 



'7(kink/sausage) — ^fkink/sausage) •5)-^'^('S), (24) 

with the superscript 'L' denoting a Laplace transform. The one-sided Fourier transform 
7] {k, Lu) of fj {k,t > 0) can now be found from equation (24), by analytically continuing s to 
icu, so that 

'7(kink/sausage) ^) = ^^ink/sausage) (k, ^<^) = G^fkink/sausage) (k, iU})F^{iu) . (25) 
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With the help of equations (19) and (25), equation (18) yields 



Akink/sausage)(A:,w)G'fkink/sausage)(^>^^)^^(^^) = ^5-(ki„k/sausage) (-A:)F^(i^), (26) 

whence we obtain 

G^fkink/sausage)!^")^^) = ^Cr(kmk/sausage) (-^)/-D(kmk/sausage) (^, t^) • (27) 

Equation (27) helps us to find G^{k,s) for any complex value of s, by means of analytic 
continuation. Thus 

^fkmk/sausage)(^' ~ ^(^(kink/sausage) (^^OZ-^Ckink/sausage) (^) ~'^"S), (28) 

which, combined with equation (24), ultimately gives the expression for the Laplace trans- 
form ^fki„k/sausage)(^'5) ^f the displacements f)(kink/sausage) (^, ^) for any complex value of s. 
This expression is given as 

^fkink/sausage)(^>s) = ^^(kink/sausage) (-A:)-F^(s)/-D(kink/sausage) (A:, -is), (29) 

from which we can find the time evolution of the displacements r/(kink/sausage)(^; ^) by means 
of Bromwich 's integral formula. Thus, 

?7(kink/sausage) {k, t) = sausage / ^^^^ <is/-D (kink/sausage) (k, -is) , (30) 

ZTTl J T— ioo 

where, the real constant T is so chosen that the singularities of the integrand lie to the left 
of the line s = T, the singularities themselves being of the nature of simple poles, that are 
the solutions of the dispersion relations £>(kink/sausage)(^, ^t-) = 0. In the present case, these 
singularities are s± — in±, with n± being simply the frequencies given in equation (20) 
above. 

To evaluate the integral given in equation (30), wc need now to consider some particular 
form of the function F{t) for the time dependent part of the external stress cr (kink/sausage) {—k, t); 
see, equation (22). A common model to choose is an exponential one, so that 
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F{t) = Aexp{-Xt), 



(31) 



with A being a constant, and the exponent A being positive definite. For this particular 
functional form of the external stress a{—k,t), equation (30) yields 



^(kink/sausage) 

— ^5'(kink/sausage) (~^)^ 

exp (s^t) 



^^■(kink/sausage) ( k) r'^+o 



exp (st) ds 



/l+C 
T-oo-D(kink/sausage)(^, S + X 

Xt) exp (s+t) 



+ 



exp I 

-'-^(kink/sausage) 

{k,-X) 

1 



+ 



^-^{kink/sausage) 



^-^(kink/sausage) (^5 



ds 



s_ + A 



S+ + X 
(32) 



where, the second expression follows from the residue theorem. The first term in this expres- 
sion is a transient term, that decays with the decay of the extrenal stress, whereas, the last 
two terms are due to disturbances that live on even after the withdrawal of the external force 
and, therefore, pertain to the eigenmode oscillations of the slab. For the sake of illustrations, 
we may consider a delta function type external stressing, so that A — >• oo, and A/A — >• 1. 
The displacements fj{k, t) then evolve as 



^(kink/sausage) (^) ^) — ^ ^^(kink/sausage) ( ^) 



exp {in^t) 



(kink /sausage) 



(fc,n) 



dn 

exp {iri-t) 



n+ 



1 - 


f^-^(kink/sausage) (^)^) 






dn 





i/C(7(kink/sausagc) 



[exp (m+i) — exp (m_i)] , 



(33) 



givng us the familiar result, that the application of an instantaneous external stress a{—k)S{t) 
on the MHD slab creates long hved excitations pertaining to the normal surface modes of the 
slab, that have frequencies n+ and n_ as given by the dispersion relations (20). Equation 
(33) further shows that the two surface modes, denoted here by a '+' and a '-' sign, are 
out of phase with the external stress a{—k,t) by 7r/2 and 37r/2, respectively. The modes 
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are of same amplitude, ie., |?7jink/„ge)(^' ^)l = l%ink/sausage)(^'^)l' but their different phase 
relations with respect to a{—k,t) lead to a difference in their respective energy absorption 
rates, as will be presented in the next section. 

3.3. Energy Absorbed by the Modes: Equation (33) shows that the rate of absorption 

of energy per unit area by each of the ^^ink/sausage) ^) uiode, from the external stress 
c^(kink/sausage)(-A:,t) is givcu by 



^ ('7(kink/sausage) (^5^)5 ^(kink/sausage '^(kink/sausage) ( ^) ^)^(kink/sausage) 

{k,t) 

exp{(m±-A)0, (34) 



kA \o'(]<imk/sa.usa,ge)i^ 




n± 




"9£'(kink/sausage) ' 




in± + A 




dn 


n± 



whence we calculate the total energy (per unit area) absorbed by the mode, by integrating 

equation (34) over a time t that is much longer than the decay time of the stress, ie. t 3> 1/A. 
This absorption is given by 



A/i(kink/sausage) (^i) = / h^^^ ('7(kink/sausagc) ik,t), (kink/sausage) ^)) dt 

J 

^|5'(kink/sausage) (~^) P''^± A'^ 



(35) 

dn \n± 

which, after substitution in favour of the modal amplitude ^(kiik/sausage)(^' ^) equation (32), 
ultimately yields 

*, /X / <9-D(kmk/sausage) (^) ^) \ 1, (±) /, ,x|2 /qc\ 

A/l(kink/sausage)(n±) = In ^ 1 T N(kink/sausage) ^) I • (36) 

We here note that, inspite of the complexities presented by the magnetic field, the ex- 
pression of wave energy presented in equation (36) is essentially the same as the one given 
in Cairns (1979) in the case of a purely hydrodynamic system. Unlike Cairns (1979), whose 
method was intuitive (see Section 1 for details) , we here derive our results directly from the 
equation of motion (18) of the system. The particular models of external stress that we 
assume for the purpose of demonstrations, do not have any bearing on our final result in 
equation (36), thus implying that this expression for the wave energy is truely a generalized 
expression for any hydrodynamic or magnetohydrodynamic system. 

15 



3.4. Waves of Negative Energy (NEW): Consider, for simplicity's sake, a frame of reference 
co-moving with the external medium, so that Woe — 0. Consider further, that the velocity of 
the slab Uqo is increased gradually from zero through positive values in this frame of reference. 
Equation (20) in Section 3.1 then shows that, as long as itoo < Ck// (aof)^^^, both the '+' 
and the modes present oscillatory surface waves with n+ < and n„ > 0, thus implying 
that the '+' wave propagates in the positive X-direction, whereas, the '-' wave propagates 
in the negative X-direction. As the value of uqo is increased through Ck// (cto/)^''^, so that 
Ck// (aof)^^^ < uqo < c^jl {aofO.efY^'^ 1 oscillatory surface modes still pertain, but now with 
n± < 0, thus implying that both the and the '-' waves now propagate in the positive 
X-direction. Thus, with the increase of the slab speed past its critical value c^// (ao/)"*^^^! 
the '-' surface wave reverses the direction of its phase propagation to be simply carried by 
the flow. In other words, the '-' surface wave changes its character from a forward wave to 
a backward wave (eg. Ostrovskii et al. 1986). To examine the energetics of these surface 
waves, we note from equation (36) that. 



A/i(kink/sausage)(n±) = ±n± (n+ - n_) - (/c, t) ^ (37) 

so that. 



A /i (kink/sausage) (^±) > 0, wheu Uqo < Ck// (^o/) , (38 a) 

and 

A/l(kink/sausage) (^+) > and A/i(kink/sausage) ('^-) < 0, 

when Ck// (ao/fte/)^^^ > Uq^ > Ck// {oioff^ ■ (38 b) 

The backward '-' surface wave is then also a negative energy wave in this particular reference 
frame - a result, that is in agreement with Cairns (1979), Ostrovskii et al. (1986) and Ryutova 

1 /2 

(1988). As the velocity rtoo of the slab passes through its threshold value Ck// (ao/ttc/) 
for K-H instability^ an unsatble region is produced by a coalescence of the positive and the 
negative energy modes. We may note that, although the sign of energy of the modes depends 
on the choice of the co-ordinate frame, the existence criterion of the unstable branch (see 
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equation (21)) is independent of such a choice. Also invariant is the total energy absorbed by 
the system from the external perturbation. For a (^-function type perturbation, this energy 
is given by (see equation 35) 

^^(kink/sausage) — ^^(kink/sausage) (^+) + ^^(kink/sausage) (^—) — ^|5'(kink/sausage) (~^) I ; (39) 

which is less than the excitation energy A/i(kink/sausage) (^+) of the '+' mode alone in the above 
example. This extra energy is released during the process of excitaion of a negative energy 
wave, thus exciting simultaneously a positive energy wave through the mode interactions in 
the presence of the external stress a{—k,t). 

4. Effects of viscosity on the v^^avemodes 

To examine the effects of viscous dissipation on the surface modes of the slab, we first observe 
that a canonical form of the stress-free equation of motion of the slab can be obtained by 
substituting d/dt for in in equation (20). This equation is 

V (kink/sausage) 

{k, t) + 2ikU j7(kink/sausage) ^) + ^(kink/sausage) ^) = 0> (40 a) 

with 

U = {a^fUoe + tto/Woo) , and {Snf = elf - a^fa^f (woe - ^^oo)^ , (40 b) 

in a frame 's'(say), in which the two fluids 'o' and 'e' are seen to move with velocities 
Uqq and Moe, respectively. In this frame of reference 's', a wave profile ri^{x,z',t) at any 
point z = z' inside the slab is seen to have a dependence given by ri^{x,z',t) = r]^{t = 
0) exp |i (rf_^t + in which the frequencies rf_^ have a nett drift term klJ, so that rf^_ ^ 

—rf_, signifying that the reflection symmetry is lost. 

Consider now a frame of refernce 'r', that moves with a relative velocity U with respect 
to the 's' frame, so that the transformation r— >s is given by 

^(kink/sausage) (^' ^) — ^(kink/sausage) (^' ^) ^^P (~'^^^^) ' (41) 
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with ?7(kink/sausage)(^' ^) Satisfying the equation of a simple harmonic oscillator 



^(kink/sausage) ^) + ^^(^'^)^^(kink/sausage) ^) " 0" (^2) 

A wave profile has a dependence given by ri^{x, z', t) = iq^it = 0) exp |i {jf-^t + kx^ | in this 
reference frame 'r', which yields n!,. = — nL — Sn (see equations (20) and (40), with U — 0). 
The slab waves in this 'r' frame thus possess a reflection symmetry, since Sn is an invariant 
that depends only on the relative velocity \uQe — uqo\, and not on the drift velocity U. 

While examining the effect of viscosity on the surface modes of the slab, we must begin 
our investigations by calculating the viscous dissipation as is seen in the reference frame 
'r'. This approach is in agreement with Rayleigh (1883), who argued that in a moving 
stream flowing with a velocity U, the nett pressure fluctuation due to viscous drag must be 

— Xd(4> — Ux^, as measured in the frame of reflection symmetry of the perturbations, 
with Xd being a viscous drag coeffecient and (l){x, z, t) being the velocity potential. In the 
present case of the magnetic modes of the slab, this requires that the rate of mechanical 
energy dissipation of the slab due to viscous damping (eg. Landau and Lifshitz 1959a) is 
given by 



^!dz+ j dz 

CO J d 



dx, 



where. 



* = 2poet'e 



(43 a) 



(43 b) 



with Uc being the kinematic viscosity in medium 'e', whereas, medium 'o' is taken to be 
inviscid. It is possible to use the classical gas dynamical formula for viscous dissipation (as 
given in equation (43)), while retaining the velocity discontinuities at the slab interfaces, 
only in such a situation, where either the internal or the external medium alone has viscous 
dissipation and the other medium is inviscid, ie., the details arising due to boundary layer 
may be ignored and the tangential discontinuity of the velocities at the interfaces still remains 
a valid condition. We however note that, ideally one should consider an anisotropic viscous 
stress tensor in the presence of a magnetic fleld (cf. Braginskii 1965), as in the situation 
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considered here. Whatever the case may be, the specific choice of the viscous stress tensor 
is not expected to change the overall stability properties of the modes about which we are 
mainly concerned in this paper. 

We now substitute the expressions 



Sul^{x, z, t) 

and 

for the various perturbations in equation (43). 
see Section 2.1), we thus obtain 



(44 a) 
(44 b) 

Using the definitions given in equation (3, 



= ^kfj (k, ^d, t) exp {±k (z ± d)} exp (ikx) 
— =F (^^) v'^ik, ^d, t) exp {±k {z ± d)} exp (ikx) , 

= ff{k,^d,t) exp{± {z ± d)} ex.p{ikx) 
— (5n) ffik, ^d, t) exp {±k (z ± d)} exp (ikx), 



(^(kink/sausage) (^)^))''7(kink/sausage)(^'^)) — —'2ipOe'^ek{Sn) ?7(kink) (^5 ^)^(kmk) (^) 

+^(sausage)(^'^)^(sausage)(^'^) ' (^5) 

for the rate of viscous dissipation in frame 'r'. 

The energy thus dissipated in frame 'r' gives rise to an increament 6Soe in the entropy 
of the system, that must be invariant in all frames. Noting that, Sn in equation (45) is an 
invariant (see, equation (40b)), we obtain the rate of increase of entropy (or, the heating 
rate) in terms of the quantities defined in the frame 's'. Thus, substituting Sn — n + kU, we 
have 



s s 

^Oe-^ (^(kink/sausage) (^' ^) ' ''/(kink/sausage) (^' ^)) = 2ipoeiyek{n + kU) [?7(kink)(^5 ^)^(kink)(^5 ^) 

+^(sausage)(^)^)^ (sausage) (^>^) > (^6) 



where, Toe is the equlibrium temperature of the medium 'e'. 

The heating rate being thus known, we demand that, the thermodynamic potential $o {— 
H—ToeSoe) must be minimum at all instants for the wave propagation to be a manifestation 
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of the system's response to its departure from equilibrium (cf. Landau and Lifshitz 1959b; 
Glansdorff and Prigogine 1971). This shows, 



S%{t) = ^oit + St) - Mt) 
dH 



(kink/sausage) 



"^(kink/sausage) 



dH 



SSi 



^(kink/sausage) ^Oe" f^— } St — 0, 



9t] (kink/ sausage) 



St 



(47) 



for any infinitesimal St. In equation (47), we have dropped the superscript 's', still indicating 
the observer's frame. With the help of equation (18), we then obtain (after substituting 

5"(kink/sausage) ( k,t) 0), 



2sinh {2kd)±4:{kd) 



8 cosh^ (kd) 



sausage ) (^5^) ^(kink/sausage) 

{k,t) 



+ 



POe (n + kUQeY 



3k'B, 



2d2 
Oe 



-Poe^'ek'^ 



16n 

2sinh {2kd)±4:{kd) 
8 cosh^ (kd) 



+ 



Poo ("- + ^^too) - 



2 k'B? 



Oo 



47r 



fik) 



?7(kink/sausage) {k, ^)?7(kink/sausage) ^) 



' I ^(kink/sausage) V/^' '^J — 



n 



(48) 



which is true for all ?7(kink/sausage)(^' Substituting d/dt — in, and also requiring a non- 
trivial solution, we then obtain the dispersion relation 



POe 



{n + kuoe) - k CAe + Poo {n + kuoo) - k cao f{k) = ipoeJ^ek {n + kU), (49) 



for the viscous surface modes of the slab. Note the factor (n + kU) in the damping term of 
equation (49). This factor differentiates equation (49) from the earlier results (cf. Kikina 
1967; Weissman 1970; Cairns 1979; Ezerskii et al. 1981; Ostrovskii and Stepanyants 1982; 
Ostrovskii et al. 1986; Ruderman and Goossens 1995), in which the viscous drag force was 
proportional to the frequency n of the waves in the observer's frame 's', thus depending on 
the velocity U of the material. In view of this important difference, it is here pertinent, that 
we discuss the significance of equation (49) in some detail. 
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The flow of tlie two fluids creates a momentum flux pooWoo + PocUoe per unit volume in 
the observer's ('s') frame, that is equivalent to imposing a velocity U — (cco/iioo + c^e/'^oe) on 
all matter in the wave profile. Moving to any other frame, where velocities of the fluids are 
Uq^ — Uqo + u and Uq^ — Uq^ + u, we have U' — ao/^oo + cte/^toe — U + u. This nett velocity of 
the wave profile appears purely due to Galilean transformation, and should not contribute 
to any process of exchange of energy or momentum within the system, and thus cannot 
contribute to dissipation. Contrary to the earlier results, the expression for the viscous drag 
must, therefore, have no explicit dependence on the drift velocity U, as is evident by the 
appearance of the invariant factor 5n — n + kU in the damping force in equation (49). 

Returning to the modes of oscillations of the magnetic slab, equation (49) yields solutions 



U± 



C^f — tte/tto/ {Uoe — ^Oo)' 



2 2; 2 



1/2 



(50) 



thus showing that, for flow velocities below the threshold for the K-H instability, the principal 
effect of viscous dissipation is to introduce a damping for both the positive and the negative 
energy modes of the magnetic slab - a result, that is in contradiction to the earlier results 
(see the references above), which predict a dissipative instability for the negative energy 
waves. The main consequence of our considering the correct Galilean transformation, while 
examining the viscous effects on the slab waves, is then the establishment of the fact that, the 
stability property of the modes remains non-singular in the presence of a small dissipation, 
that changes only slightly the threshold for the K-H instability of the slab, with the modified 
instability criterion given by 



koe - I^Ool > 



1 a 



1/2 



(51) 



that smoothly approaches the adiabatic criterion in equation (21) for a vanishingly small 
kinematic viscosity Ue ^ 0. 
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5. Concluding remarks 



Occurrence of magnetic structures is abundant in various astrophysical situations, such as in 
the solar photospheric flux tubes, or in the solar coronal plasma loops. Such magnetic struc- 
tures are often associated with field-aligned plasma flows, with the velocities of these flows 
being different inside the structures than those outside, thus producing shearing motions in 
the plasma medium. Detailed understandings of the complex interactions of such shearing 
flows with the oscillatory motions of the structures are necessary to study accurately the 
energy transport processes in astrophysics, such as the mechanisms of non-thermal energy 
transport from the solar sub-surface layers to the upper solar atmosphere. Certain investiga- 
tions have been carried out (eg. Ryutova 1988; Nakariakov and Roberts 1995; Nakariakov et 
al. 1996; Ruderman and Goossens 1995; Ruderman et al. 1996; Joarder et al. 1997) in this 
direction, that highlited the role of negative energy waves in such processes. As a further 
contribution to such investigations, we here examined in detail certain specific aspects of 
the interactions of magnetohydrodynamic waves with shearing flows, and particularly of the 
negative energy waves, by using a self-consistent thermodynamic approach. This approach 
helped us to generalize the expression for the hydrodynamic wave energy given in Cairns 
(1979) to magnetohydro dynamics (see, equation (36) in Section 3.3), thus enabling us to 
calculate the energy of the hydromagnetic waves (of course in the harmonic approximation) , 
when the linear dispersion relations of such waves are known along with the observationally 
obtained informations regarding the wave amplitudes. Once the wave energy is thus cal- 
culated, equation (35) then guides us to obtain a rough estimate of the generating stresses 
a{—k,t) of the waves. Such estimates of the stresses may be of great importance in sev- 
eral astrophysical situations, particularly in solar MHD cases, where such estimates may 
provide us with some clues regrading the physical processes that may be taking place in 
the sub-surface layers of the Sun, or in the regions of complex magnetic topology in the 
solar atmosphere, about which we have very little direct observational evidence. Finally, by 
incorporating viscosity, we obtain the dispersion relations (equation 49 in Section 4) which, 
while precluding the possibility of dissipative instability,- sets the correct conditions for the 
stability of the system (equation (51)) and also yields the time scales for the decay of the 
disturbances in the surface modes of MHD systems. It is to be hoped, that the present study 
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would provide us with some guidance in gaining further physical insights into the complex 
nature of the interactions between the magnetic field and the fluid flows in various astro- 
physical systems,- both for estimates in terms of energetics as also in the study of evolutions 
of MHD eigenmodes. 
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Appendix. Energy density of an incompressible MHD plasma. 



The energy density of an incompressible magnetohydrodynamic system is defined as 

e{r,t) = s,{r,t) + ?^^ + ^u'{r,t), (A.l) 

where, is the thermodynamic internal energy of the plasma, and the quantities po, u and 
B are as defined in equation (2) of the main text. 

In the present case, we consider the fluid to be incompressible, and the hydrodynamic 
processes to be adiabatic. In that case, we write (cf. Landau and Lifshitz, 1959b) 



£i(r,i) = fi{r,t)/v, (A.2) 

where, /i (r, t) is the local chemical potential of the system, and v is its specific volume. Thus, 
for any fluctuation in the thermodynamic state of the system, we have 



Ssi (r, t) = 6fi (r, t) /v, (A.3) 
where, v is a constant in an incompressible fluid. Further, 



Sfj, = -sdT + vdp, (A.4) 

where, s is the specific entropy of the system. 

If we now assume the electrical and the thermal conductivities of the fluid to be infinite, 
and its viscosity coefficient to be zero, then the system cannot support any thermal gradients, 
and, therefore, ST is zero at all points. Combining equations (A.2-A.4), we then find that 
the fluctuations in the thermodynamic energy density to be (see, Sommerfeld 1950 for an 
alternative interpretation) 



Sp, 
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(A.5) 



so that, neglecting the integration constant, we ultimately obtain equation (6, Section 2.2) 
of the main text, ie., 

e (r, t)^p (r, t) + + (r, t) (A.6) 

for an incompressible fluid. 
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